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fS) . the other galaxies of the cluster. This influence manifests itself only as tides in standard 



Newtonian gravity, meaning that the systematic acceleration of the centre of mass of the 
galaxy has no consequence. However, in the context of MOND, a consequence of the 



own self-gravity of the galaxy. We show here that, in this framework, initially axisymmetric 
elliptical galaxies become lopsided along the external field's direction, and that the centroid 
of the galaxy, defined by the outer density contours, is shifted by a few hundreds parsecs 
with respect to the densest point. 
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jrt ' Abstract: We point out an interesting theoretical prediction for elliptical galaxies residing 

inside galaxy clusters in the framework of modified Newtonian dynamics (MOND), that 
could be used to test this paradigm. Apart from the central brightest cluster galaxy, C ^ 

other galaxies close enough to the centre experience a strong gravitational infiuence from ■>«—-' 



o 
o 



C^ ' breaking of the strong equivalence principle is that the systematic acceleration changes the v — J 
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1. Introduction CD 

Modified Newtonian Dynamics ||l|, Q] is a recipe linking tlie gravitational field in galaxies 
with their baryonic content. This recipe makes many successful predictions on galaxy 
dynamics, such as the baryonic Tully-Fisher relation [y], and the precise shape of rotation L^-' 

curves of both high and low surface-brightness galaxies in the field [Q], as well as those of C !) 

tidal dwarf galaxies (e.g., Gentile et al. 2007). It could also hold clues to the observed 
universality of dark and baryonic gravities at the effective dark halo core radius |5|, y, Q- 
The successes of this recipe in galaxies could be an emergent phenomenon, linked with the 
complex feedback between baryons and dark matter in the process of galaxy formation Q , 

but a more radical explanation of these successes is a modification of gravity on galaxy C J 

scales. This interpretation requires to postulate additional dark matter at extragalactic C ^ 

scales (e.g., Sanders 2007; Angus et al. 2007; Angus 2009), but let us note that most ^~^ 

relativistic versions of such modified gravity theories inevitably include dark fields that 
seed structure growth, and may play the role of dark matter in galaxy clusters [P]. This 
reduces calculability on large scales, and means that the modified gravity interpretation of 
the MOND recipe is more easily falsifiable on galaxy scales. Falsifying such an alternative 
approach to the problems of galactic dynamics is actually needed in order to more firmly 
anchor the current cold dark matter (CDM) paradigm as the one that should ultimately 
be confronted with data at the sub-galactic, galactic and cosmological scales. 

One of the unique properties of the modified gravity version of MOND in galaxies, 
compared to its Newtonian CDM counterpart, is that it violates the Strong Equivalence 
Principle when considering the external gravitational field in which a system is embedded. 
Hence, the effects of the gravitational environment are important, and fundamental con- 
cepts such as the escape speed from a system are fully dependent on the external field. This 



can lead to extremely interesting and testable predictions. For instance, [10| analytically 
predicted in the Solar System an anomaly due to the external field of the Galaxy, in the 
form of a quadrupole field that may be detected in future measurements of planetary and 
spacecraft motions. On the other hand, |11] and [12, ^] studied the effect of the external 



o 



field from large scale structure on the Milky Way and showed that the MOND predictions 
for the galactic escape speed were quite realistic, while allowing the Large Magellanic Cloud 
to be on a bound orbit. 

Galaxies inside clusters present a very interesting environment where the external field 
has a significant influence. These galaxies experience both their self-gravity and the gravity 
of the other galaxies of the cluster (except for the galaxy at the centre of the cluster): 
the latter manifests itself only as tides in standard Newtonian gravity, meaning that the 
systematic acceleration of the centre of mass of the galaxy has no consequence. But in the 
context of MOND, a consequence of the breaking of the strong equivalence principle is that 
the systematic acceleration changes the own self-gravity of the galaxy. 

In this Letter, we emphasize that, in MOND, the combination of the self-gravity with 
the external acceleration breaks the front-back symmetry of a purely self-gravitational 
system. As elliptical galaxies are abundant in clusters, they provide an ideal case to test 
this front-back asymmetry. Offsets such as those discussed in this letter have also been ' \J 

studied in the recent literature in the context of MOND N-body simulations |14], CDM (^ ^ 



on the left-hand side of the modified Poisson equation [^, 12, 13 1: because the external 
field has a direction, this total gravity is increased by the external field on one side of the 
galaxy centre along the x-axis, and decreased by it on the other side. This effect is maximal 
when the internal and external gravities are of the same order of magnitude. Clearly, for an 
originally spherically symmetric distribution, the potential and its derivatives are then still 
axisymmetric, but the same conclusion holds when applied to an axisymmetric system, and 
if the direction of the external field does not coincide with the original axis of symmetry, 
the axial symmetry is broken. We also computed a flattened (axis ratio 0.7) Plummer 
model with a diagonally pointing external field. As expected, we see the twisted and offset 
potential contours (Fig. 2). 

An additional interesting effect is the following: because of this lopsidedness of the 
potential in the outer parts, an additional "external" field is generated, which acts on the 
inner parts of the galaxy, akin to the effect in the inner Solar System pointed out by |10|| . 
This will cause the centre of the galaxy to be slightly shifted with respect to the point 






N-body simulations |15], and lensing observations [^]. ^ — v 

2. Lopsidedness and offset of the centre of gravity . 

Let us start with a simple numerical experiment: we place a Plummer sphere into an C !) 
homogeneous external field typical of a cluster environment, of the order of ~ 10^^ cm s~^, 
and we solve the modified Poisson equation of MOND, which is a non-linear elliptic partial 
differential equation Q: the boundary conditions are taken as in ||T^, and the boundary 
value problem is solved by using the iterative Newton method in spherical coordinates [|l^ . 

The isodensity and isopotential contours in the free- falling frame are then plotted on Fig. 1, C J 

where the external field is applied in the negative x-direction: clearly, one sees that the C ^ 

isopotential contours are flattened and lopsided, in the sense that they do not exactly /- — ^ 
correspond to the same radii on the positive and negative x-axis. The explanation for this 
is that, in MOND, the iota/ gravity (internal+external) enters the dielectric-like ^-function 
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Figure 2: Same as Figure 1 for an axisym- 
metric Plummer spliere with, an axis ratio of 
1 : 1 :0.7. 



Figure 1: Isodensity (dotted) and MOND 
isopotential (solid) contours of a Plummer 
sphere with a 5 x 10^" M© mass and a Plum- 
mer scale-length of 1 kpc, embedded in an 
external field of 1.2 x 10~^cms~^ with the 
direction indicated by the arrow. 



where the internal gravity is zero (see also |]14[). In the Plummer sphere model of Fig. 1, 
the internal gravity is zero at x = 20 pc. Although this is small, it of course depends on 
the model parameters, and it means that the initial configuration can never stable and 
that the photometric centre of the galaxy should shift towards this new centre of gravity, 
thereby increasing the lopsidedness of the galaxy. We hereafter investigate whether this is 
the case with numerical N-body models of realistic elliptical galaxies. Starting from initial 
conditions built with Schwarzschild modelling, we let the galaxy evolve towards its final 
state, and check whether it becomes significantly lopsided. 

3. Schwarzschild and N-body models 

As a first step, in order to get our initial conditions, we build static equilibrium models 
of axisymmetric ellipticals in MOND, with an external field of 1.2 x lO^^cms^^. The full 
details of our procedure are given in a companion paper (Wu et al. 2010, submitted). We 
use the so-called Schwarzschild method 111, 19, 20] to reproduce the system self-consistently 



^ by a linear superposition of orbits in the computed MOND potential of an axisymmetric 
|21] profile with a total (baryonic) mass M = 5 x lO^'' Mq, and an axis ratio of 0.7. Even 
when using 27360 orbital building blocks, we still find a significant residue of 4% between 
the Schwarzschild model and the input density: this means that a fully self-consistent model 
cannot be found for an axisymmetric density model embedded in a strong external field 



^Historically self-consistent galaxy models in Newtonian dynamics mean equilibrium models purely in 
the gravity of the baryons, without Dark Matter. In MOND this obviously is also the case since no dark 
matter is supposed to be present. However, the original density reproducing the external field here is not 
modeled, hence our models are not purely under self-gravity of the model galaxy. 
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Figure 3: Color maps (centroid details on Figure 4: Same as Fig. 3 after t=90 Tgimu- 
the scale of 1 kpc) of xz-projected density for The blue circle in the centre shows the ori- 
our initial conditions gin point (i.e. centre of mass), roughly corre- 

sponding to the centroid defined by the outer 
surface density contours. It is clearly offset 
from the central density distribution. The 
arrows indicate the direction of the external 
field. 

in MOND. We thus expect the system to evolve towards a non-axisymmetric equilibrium 
configuration. 

We then turn this initial Schwarzschild model into a live N-body system by Monte 
Carlo sampling of the phases of the orbits. Starting with these initial conditions, we use 



the iV-MOZ)F code |22| and let the system evolve. Figs. 3, 4, 5 and 6 show the evolution 
of projected density of the model after 0, and 90 simulation time (1 simulation time is 1 
Keplerian time at the typical scale of 1 kpc, see Eq. (6) in |23], ITgimu = 4:.7 Myr). One 
can see clearly an offset of ~ 200 pc of the central densest point from the centroids defined 
by the outer projected density contours. The combination of orbital anisotropy and the 
broken axisymmetry makes the evolution quite complex. We will report these elsewhere. 
By and large this live N-body simulation confirms our expectation that the axisymmetry 
is broken in the presence of an external field, which shows up as a twist of the isophotes 
and an offset of the density peaks. 
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4. Conclusion 

We estimate that for typical parameters of elliptical galaxies residing in large clusters, a 
centroid shift of a few hundreds of parsecs is expected in MOND, as compared to the 
centroid of the outer isodensity contours of the galaxy. While tides can also cause some 
lopsidedness in Newtonian and MONDian gravity, there is a major distinction with the pure 
MOND effect under scrutiny here. The tidal effects are small inside the tidal radius, and 
are traditionally neglected when dealing with dense elliptical centres. The fact that these 
effects can be neglected in our analysis is confirmed by the fact that the average density of 
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Figure 5: Same as Fig. 3 on the scale of 
10 kpc 
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Figure 6: Same as Fig. 4 on the scale of 
10 kpc. This figure shows clearly the twist of 
isophotes between the inner and outer den- 
sity distribution. The arrow indicates the di- 
rection of the external field. 
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our model galaxy inside the saddle point where the external field starts to dominate is larger 
than the cluster density. The saddle point is estimated to be 15 kpc from the Plummer 
sphere (total mass 10^^ Mq) centre (Fig.l of (ij]). Here the average density enclosed is 
^ 4x 10~^Mqpc~^, typically more than one order of magnitudes larger the typical cluster 
density at that same radius (1000 x pcrit ~ lO~'^M0pc~^). Note that [24] argues that such 
tidal criteria holds in MOND as well in Newton. 
Thus, we conclude that: 

• If real elliptical galaxies in rich clusters show perfectly symmetric light with no sig- 
nificant offsets between the centroids of the inner and outer contours, the classical 
version of MOND is likely excluded. This predicted lopsidedness of galaxies inside 
distant rich clusters should be falsifiable with photometry only, using, e.g. the VLTI. 
To resolve a centroid offset of 200pc in an elliptical galaxy in a rich cluster of typical 
internal gravity ~ 10~^ cms~^ at a distance of 160-210 Mpc (e.g., Abell 1983, Abell 
2717, MKW9, [P5| , Pq ] would require a minimum angular resolution of 0.5 arcsec. The 
offsets at the centroids of galaxies are between 10 pc and 500 pc: for an elliptical in 
Coma (100 Mpc far away) this corresponds to 0.02 to 1 arcsec, while for an elliptical 
in Virgo (20 Mpc far away) it corresponds to 0.1 to 5 arcsec. 

• On the other hand, observing such a lopsidedness would not consitute a direct fal- 
sification of the CDM paradigm. Indeed, CDM simulations predict that isolated 
CDM halos are typically lopsided due to a lack of relaxation, so this should be the 



case inside clusters too. [27] studied the shapes of a large sample of ten thousand 
pure CDM halos, and found that any offset between the barycentre (CoM) and the 
density-centre (potential minimum) is correlated with a massive satellite inside an 
unrelaxed halo, i.e., the offset is a measure of unrelaxedness in the CDM context. 
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This is characterized by an offest of order of 1 or 2 percent between the central dark 
matter density maximum and the centre of mass at the virial radius. Let us note 
that if this can be translated into a predicted lopsidedness for the light distribution 
of isolated ellipticals in CDM, this would be a clean test since MOND does not pre- 
dict such an offset (except at extremely large radii) for isolated galaxies. What is 
more the centres of some ellipticals are home to strong star formation which could 
also distort the observed central isophotes. The isophotal distortions predicted in 
MOND, if actually observed, could thus be due to this. However, for galaxies inside 
clusters, these more mundane offsets would not be expected to be aligned with the 
direction of the local external field pointing to the cluster centre. If a statistically /- 

relevant preferred direction is found for the lopsidedness of galaxies residing in clus- | , 

ters, it would thus be in support of MOND, although it would not constitute a direct j~^ 

falsification of CDM. M 

Nevertheless, a caveat is that our simple models here were based on a constant external /- — n 

field of ~ 10~® cm s^^. However, in reality, while galaxies are orbiting inside clusters, the 
external field acting on them varies with time, both in direction and amplitude. So the 
next step will be to make the external field vary as a function of time in our simulations, '""^ 

to check how this affects the predicted lopsidedness. L\J 

Finally, we also note that the flattened and lopsided potential created by the external ^ ^ 

field would also be valid in disk galaxies, although most of them are in the field or at the 
outskirts of galaxy clusters. Still, the flattened MOND potential would create a differential 
force with a component normal to the disk, hence a specific torque. This could cause dif- 
ferential precession of the disk angular momentum vector, and could lead to the formation 
of warps, even in a very low external field for isolated galaxies |28|. C J 
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